Materials design from non-equilibrium steady states: driven graphene 
as a tunable semiconductor with topological properties 
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Controlling the properties of materials by driving them out of equilibrium is an exciting prospect 
that has only recently begun to be explored. In this paper we give a striking theoretical example 
of such materials design: a tunable gap in monolayer graphene is generated by exciting a particular 
optical phonon. We show that the system reaches a steady state whose transport properties are the 
same as if the system had a static electronic gap, controllable by the driving amplitude. Moreover, 
the steady state displays topological phenomena: there are chiral edge currents, which circulate a 
fractional charge e/2 per rotation cycle, with frequency set by the optical phonon frequency. 



Non-equilibrium quantum systems constitute a natu- 
ral frontier in physics that is only beginning to be probed 
by theory and experiment. Non-equilibrium methods can 
be used to study [T] and control [2 ^ the properties of 
condensed matter systems. Particularly exciting is the 
possibility of engineering the properties of novel materi- 
als, like graphene, by driving them out of equilibrium, 
paving the way for applications to devices. Graphene's 
gaplessness poses a critical challenge to such applications, 
as the development of graphene-based semiconductors is 
predicated on the ability to induce a gap. 

In this Letter we present a theoretical study in 
graphene of this non-equilibrium approach to materi- 
als design. We demonstrate the possibility of inducing 
a gap in monolayer graphene by the excitation of op- 
tical phonon modes. The gap is controlled by a time- 
dependent Kekule-pattern bond density wave, which ap- 
pears in the effective field theory as a complex-valued 
order parameter A that rotates with the frequency fl of 
the driven phonon mode. The time-dependence in this 
order parameter is completely removable by an axial (val- 
ley) gauge transformation, which can be viewed as a kind 
of "boost" to a co-moving "reference frame." The gauge 
transformation has no effect on the coupling of the system 
to a heat bath, thereby guaranteeing thermal equilibra- 
tion in the new frame, and leaves the fermion currents 
invariant. This implies that the electric response of the 
system is equivalent to that of one with a static gap; all 
non-equilibrium aspects of the problem are removed and 
the system can be studied as if it were at equilibrium. 

The topological consequences of the Kekule gap have 
been studied in the static case, revealing that fraction- 
ally charged states can emerge that are bound to vor- 
tices in the order parameter A 1^. In the driven case, we 
show that further topological phenomena arise: the sys- 
tem supports chiral edge currents of magnitude Jcdgo = 



e ^I/Att, while the current in the bulk vanishes. These re- 
sults suggest the possibility that driven graphene could 
be used as a tunable semiconductor with nontrivial topo- 
logical properties. 

Let us consider spinless electrons hopping on a hon- 
eycomb lattice A according to the time- dependent tight- 
binding Hamiltonian 
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I] + ^tr,j{T)] albr+s, + h.C. , (1) 

where r is time and aj. and ^^^^ . are fermionic creation 
operators at sites r E Aa and r + Sj G As , with A^ and 
Ab the two triangular sublattices forming the hexago- 
nal lattice A. The vectors Sj (j = 1, 2, 3) connect a site 
r £ Aa to its three nearest neighbors at r + Sj E Ab 
located a distance \sj\ = d away. The uniform hopping 
amplitudes t are modulated by time- and site-dependent 
perturbations Strj{T). In the absence of such perturba- 
tions {5trj{T) = 0), the Hamiltonian ([T]) can be diago- 
nalized in momentum space, and the single particle spec- 
trum has two Dirac points at fc = K± = zfc ^^^ (li 0)- 

We shall now consider the perturbations Str.jir) that 
result from the excitation of the highest-energy optical 
phonon modes at wavevectors K± with frequency fl. The 
atomic displacements from the lattice sites va.b G Aa.b 
are 

(r-A,i3,r) = c± e^-A.-^ie-^"^ «f c.c. (2) 

The coefficients c± are the amplitudes of the excited 
waves. The normal mode vectors u^.' for the highest- 
energy optical modes with frequency at wavevectors 
K± can be determined from a classical analysis of the 
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lattice displacements [SI |S] and are given by 



2 \Ti 



and 



2 \±i 



(3) 



To determine the form of the hopping modulations 
StrjiT) resulting from the phonons, we consider the 
changes in bond lengths due to the atomic displacements 
([2]) when either the mode at or K- is excited. For 
small displacements, the change in the length drjir) of 
the bond connecting site r and r + Sj is [7] 
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K± 



u 



K± 



(r- 



(4) 



Substituting ^ and ^ into Q and using e'^± ''^ = 
g±i^(j-i) Qj-^g obtains 



d 



where the vector G — Kj^ — K- — 2K+ connects the 
two Dirac points. The modulation in the hopping am- 
plitude is related to the change in bond length through 
Str.j{T)/t = a 6d^ j{T)/d, where a ~ 3.7 is the di- 
mensionless electron-phonon coupling [7 . The resulting 
Str.j^r) can be written as 



Str,jiT) 



A(t) e^-f^+-^^e*<^"' + c.c. 



where 



A(r) = 



for the 
for the K_ 



mode 
mode. 



(6) 



(7) 



The hopping modulations Q have the form of a Kekule 
distortion with an order parameter A(t) [4 that is time- 
dependent. Therefore, exciting either the Kj^. or the 
K- mode independently yields a Kekule order param- 
eter that rotates in time with frequency f2 in opposite 
directions for the two modes. 

Without loss of generality, we henceforth consider the 
case where the mode is excited, and write A(r) — 
|A| e*'^(^), where (^(r) = Qt + ip. AU the results for 
the K mode are obtained from those below by taking 

n -n. 

We study the consequences of this rotating order pa- 
rameter in the context of the effective Dirac field the- 
ory of the system, which is valid in the limit where the 
fermions have relativistic (hyperbolic) dispersion. In or- 
der to ensure the validity of this approximation we re- 
quire |A| /t <C 1 and ^ 1, where the uniform hop- 
ping amplitude t sets the kinetic energy scale of the prob- 
lem. In this regime the Hamiltonian ([I]) corresponds, 
to first order in a gradient expansion, to the Dirac La- 
grangian density [H |H] 



C = ^ 7^(i5^+75Ap)-|A| e- 



(8) 



with /X = 0, 1, 2, 5' = vjit^o and 4x4 Dirac matrices 
.,0 _ /O 
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where 1 is the 2x2 unit matrix and ct^ are the three Pauli 
matrices. The Dirac spinor — (6p.4. ap_|_ a}p _ b^p-) 
collects the creation operators aj, _|_ and 6p _|_ for the ± 
species on sublattices A and B, respectively. The ax- 
ial gauge field ^5^, examined in a different context in 
Ref. 8j plays an important role in the discussion of the 
asymptotic steady state of the driven system. The spa- 
tial components ^5 i correspond physically to acoustic 
phonons and strain in the graphene lattice. If the lattice 
is strained uniaxially, the hopping amplitudes change, 
and the Dirac points shift away from K±. In this case, 
the A^i acquire a non-zero average value. In addition, 
acoustic phonons, either in-plane or out-of-plane, dynam- 
ically stretch the bonds, leading to fluctuations of Ac,i 
around the average. These acoustic phonons provide a 
thermal bath and their coupling to the electronic degrees 
of freedom provides a system-bath interaction, which en- 
ables the system to reach an out-of-equilibrium steady 
state. 

We now observe that the time-dependent mass term 
in the Lagrangian (|8| can be made constant by the axial 
(valley) gauge transformation 



2" 
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(9) 



where i — 1,2. The transformed Lagrangian is found to 
be 



£ = § 7^(*a^+75i5M)"|A| 



(10) 



where we used {75,7^} — 0. This transformation maps 
the problem to a frame of reference which is "co-moving" 
with the Kekule mass, so that the Lagrangian is no longer 
explicitly dependent on time. 

The vector current operator j'^ = ^'7'^^!', which is as- 
sociated with the electric response of the system, and 
the axial current operator j'g — ^^^^^^^'^ are invariant 
under Q. Furthermore, the spatial components A^i of 
the axial gauge field are also invariant under ([9]). Since 
we have taken the fluctuations in Ac, i to act as a heat 
bath, we conclude that this transformation leaves the 
bath invariant. Moreover, it also leaves the system-bath 
coupling A^ijl invariant. Therefore the transformation 
([9]) removes all time-dependences — those of the system, 
the bath, and the system-bath interactions. The remark- 
able consequence is that the non-equilibrium steady state 
of the time-dependent system corresponds to a thermal 
equilibrium state in the co-moving frame! 
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Consequently, the Hamiltonian T-L corresponding to the with < 6* < tt and < < 27r. Equations ( 14 1 and ( 15 ) 



transformed Lagrangian ( 10 1 can be analyzed in the time- 



independent Schrodinger picture at thermal equilibrium. 
% takes a particularly simple form in the absence of 

= 0, i.e. A50 = — r2/2 and 



strain, in which case 
is, = 0: 



H 
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|A| e"^ 1 
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(T ■ p 



(11) 



where cr is the 2D vector of Pauli matrices and p — — zV. 
The eigenvalue problem Hip = Eip has been solved in [9j 
in the context of the superconducting proximity effect in 
topological insulators [TU] ; the four energy eigenvalues of 
the Hamiltonian (11) are given by 



E. 



(12) 



Evidently the gauge transformation (|9| maps the time- 
dependent problem of Eq. ([s]) to a time-independent 
problem with an energy gap 2 | A | . 

It is important to observe that, because the vector 
current operator is invariant under ([o]), all observ- 
ables associated with j'' can be calculated from the static 



Lagrangian (10) without dealing with the original time- 



dependent mass. In particular, the conductivity tensor 
Gij obtained from the Kubo formula written in terms 



of the current operator can be computed from ( 10 1 



Consequently, the driven graphene system effectively be- 
haves as a semiconductor with a gap 2 |A| tunable by the 
amplitude of the optical phonon mode. 

We shall next demonstrate that the rotating Kekule 
mass in the Lagrangian ([8| gives rise to topological phe- 
nomena beyond those that have been found in the static 
case. To do this, we follow [111 in studying a variant of 



7V]*, (13) 



where the scalar field /i = iJ.{x) corresponds to a stag- 
gered chemical potential that establishes an energy im- 
balance between the sites of and A^. The Kekule 
field A = |A(a;)|e*'^(^'^), where 0(a;,T) = VLt + ip{x), 
now carries an explicit spatial dependence. The fields 
/i and A correspond to independent masses in the La- 
grangian (13), i.e. the total effective mass of the charge 
carriers is \J ^ I^P- The vector current density in the 
presence of (space- and time-dependent) masses and A 
is given by [11] 

(J^) = e^e'^"^ {do.X*dpX " ^ [(1 - 2\x\^)A,p] } , 

(14) 

where e is the electron charge, e'^"^ is the Levi-Civita 
symbol, and the complex-valued auxiliary field x ^ 
sin(6'/2) e'^, where 

cosg = . sin6> e'*^ = , „^ . . , (15) 



|A|2- 



Vm' + |ap 



form the basis of our discussion of the topological cur- 
rents resulting from the time-dependence of the Kekule 



mass term in (13 1. We use p to define an edge, setting 
/i — ^ in the bulk and using the limit \fj,\ ^ 00 to define 
an insulating region outside the sample |16j . 

The current density in (14 1 is gauge-invariant, so one 
can compute it in the reference frame where (f> has a time 
dependence or in the co-moving frame where (j) (and x) 
are time-independent. It follows that the averaged charge 
and current densities are 



{p)=e-e''^^da*d,x={p)s 
Zn 

ZTT 



(16a) 
(16b) 



where z is the unit vector perpendicular to the plane of 
the sample. 

Several observations are in order. First, the charge 
density in the case of the time-dependent Kekule mass is 
identical to that in the static case. Second, the current 
density is non-vanishing and proportional to the rotation 
frequency ft. Notice that the rotating mass breaks time- 
reversal symmetry, and therefore it is possible to have a 
non- vanishing current. Third, if |A| does not vary spa- 
tially, the current vanishes; this is the case in the bulk of 
a uniform graphene sample, where we take |A| to be con- 
stant. Fourth, there are necessarily edge currents, which 
we shall now discuss in detail. 



It follows from ( 16b I that the currents flow perpen- 
dicular to the gradient of |A|. At the boundary of the 
sample |A| must go from constant to zero. Therefore an 
edge current should flow parallel to the boundary, within 
the region where |A| varies in space, (see Fig. [1]). The 
edge current is given by 



Jc.dfra — 



(i X d£) 



(17) 



where ^ is a path that traverses the boundary. In the 
interior of the sample |A| is non- vanishing, so we can set 
^ — T' 0, and using Eq. (15) we obtain that |xinP — >■ 1/2. 
Outside the sample, |A| — > and |^| 00. Depending 
on whether > or /z < we obtain |xoutP — or 1, 
respectively. Therefore, we arrive at the edge current 

e n 



J, 



edge 



2 2n 



sgn/x . 



(18) 



The linear relation between Jodge and Q has a quantized 
coefficient. Note that because = 2tt/T, where T is 
the rotation period, the current Jedgc carries a fractional 
charge ±e/2 per rotation cycle [IT]. This chiral current at 
the boundary of the steady state bulk insulator is a topo- 
logical property of the out-of-equilibrium system; the cur- 
rents are quantized and protected against details at the 
edge, including disorder. 



FIG. 1: Chiral edge current resulting from the out-of- 
equilibrium steady-state arising from the excitation of optical 
phonons at wavevector K+. The direction of the current is 
inverted for K- phonons, for which Q. — >■ —Q.. 



The chirality of the edge currents depends on whether 
the K+ or phonon mode is excited. However, the 
chirality of the current also depends on sgn /i. We now 
offer a physical explanation of this fact. The mass ^ was 



included in the Lagrangian ( 13 ) as a means of terminat 



ing the sample with an insulating region. In a physical 
graphene flake, our findings therefore indicate that the 
sign of the edge current is determined by the specific 
shape of the sample. Notice that the direction of the 
current obtained from the field theory cannot change un- 
less IX changes sign outside the sample. But if this is the 
case, there will be domain walls separating these regions 
that support gapless modes. Indeed, these walls serve as 
quantum wires |12j attached to the sample, as shown in 
Fig. [2] The direction of the edge currents reverses at the 
contacts, as shown in the figure. Current conservation 
requires that currents of magnitude Jwiro = eQ./2TT flow 
in the wires, splitting equally at the contacts and travel- 
ing around the edges of the sample. The graphene flake 
in this scenario becomes a pump |13j that transports a 
charge e per rotation period T. 

The next observation concerns zero modes in graphene, 
which are supported in the presence of vortices in the 
order parameter A ^ . An external chiral gauge potential 
was added to render finite the vortex energies, thereby 
deconfining them [8] . Such a vortex background can also 
exist in our time-dependent scenario. In the co-moving 
frame this involves adding ^5^0 = — to the static 
problem. We find that zero-energy modes persist both 
with and without , consistent with the findings of jT4( 

Our final observation concerns the size of the gap 
that can be achieved by excitation of the optical phonon 
modes at K±. From Eq. ([t]) we obtain that |A| = 
Satjcil/d, where |c±|/d measures the relative displace- 
ment of the atoms from their equilibrium positions due 
to the phonons and is controlled by the intensity of the 
excitations. Using a ~ 3.7 and t w 2.8 eV for graphene, 
one obtains for a relative displacement |c±|/(i « 0.04% 
that 2|A| « 0.025 eV, corresponding to room tempera- 




FIG. 2: Currents in the presence of domain walls between 
regions with fi > and fi < 0. The edge currents have oppo- 
site chiralities to either side of the wires. The current pumped 
per cycle is an integer multiple of e, while a fraction e/2 goes 
around each side during the cycle. 



ture scales. 

In summary, we have illustrated a mechanism for open- 
ing a tunable Kekule gap in graphene by exciting an 
optical phonon mode at or K . This gap corre- 
sponds to a complex-valued order parameter A in the 
continuum theory that rotates in time with frequency 
n. The time dependence of A is completely removable 
by a gauge transformation which has no effect on bath 
degrees of freedom and leaves the current operators un- 
affected. The electric response of the system is therefore 
equivalent to that of one with a static gap. Furthermore, 
the system is found to support chiral quantized currents 
that are localized in regions where |A| varies spatially. 
In particular, there are edge currents whose chirality de- 
pends on the shape of the sample and on which of the 
K± phonon modes is excited. 
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